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Abstract. An upper bound for the number of Hamiltonian cycles of symmetric diagraphs is 
established first in this paper, which is tighter than the famous Mine's bound and the Bregman's 
bound. A transformation on graphs is proposed, so that counting the number of Hamiltonian 
cycles of an undirected graph can be done by counting the number of Hamiltonian cycles of its 
corresponding symmetric directed graph. In this way, an upper bound for the number of Hamiltonian 
cycles of undirected graphs is also obtained. 
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I. INTRODUCTION 

o 

Let G — (V, E) be an undirected graph with a vertex set V = {v%, v%, ■ ■ ■ ,v n } and edge set E. An edge from to 
Vj is denoted by (vi,Vj). For simplicity, the vertices are also denoted as V = {1, 2, • • • , n\. A Hamiltonian cycle of G 
is a closed path that visits each of the vertex once and only once. Similarly, if G is a directed graph, a closed directed 
path which visits each of the vertex once and only once is a Hamiltonian cycle of a directed graph. In this paper, we 
use the notation v\V2 ■ ■ ■ v n v\ and (v\, %)%,■■■ , v ni vi) to denote Hamiltonian cycles in undirected and directed graphs 
respectively. 

It is well known that the decision problem whether a graph contains a Hamiltonian is NP-complete. Hence, counting 
| the number of Hamiltonian cycles of a graph is a hard problem too. G.A.Dirac0] shows the existence of Hamiltonian 
O . cycles in the undirected graph G of minimum degree at least (1/2 4- e)n, where n is the number of vertices in G 
and e > 0. Counting the number of Hamiltonian cycles in such graphs is still #P-complete[8]. Hence algorithms 
and analysis are developed for approximating or estimating the number of Hamiltonian cycles in both directed and 
' undirected graphs . The best asymptotic result of the number of Hamiltonian cycles on random graphs is obtained 
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, by Janson 1 1 1 . N.Alon et. al. show better lower and upper bounds of the maximum number of Hamiltonian cycles 
in an n-tournament problem [l], [2j. This paper mainly focuses on bounding the number of Hamiltonian cycles. An 
upper bound of the number of Hamiltonian of arbitrary symmetric directed graph is proposed. We also prove that 
the number of Hamiltonian cycles of an undirected graph equals half of that of its corresponding symmetric direct 
graph. Hence, any bound on a directed graph can be directly applied to bound the number of Hamiltonian cycles of 
an undirected graph. 

Our novel bound on a symmetric directed graph is better than one of its natural bounds Mine's bound and tighter 
than Bregman's bound in many cases, e.g. when out-degrees are bounded by a constant K < 5. Our proof of the 
new bound is mainly based on a random algorithm for counting the number of Hamiltonian cycles on a directed 
graph, which is modified from Rassmussen's algorithm [15]. To apply the result on a symmetric directed graph to the 
undirected graph, a very simple but useful transformation that transforms counting the number of Hamiltonian cycles 
of a undirected graph to that of a symmetric directed graph is proposed. 

The structure of this paper is as follows. Some nature bounds for the number of Hamiltonian cycles led by matrix 
permanent are introduced in section II. The Rassmussen's algorithm for counting the number of Hamiltonian cycles 
is discussed and a modified algorithm is presented in section III. Some fundamental properties of the algorithms are 
given. A new bound on a symmetric directed graph is presented in section IV. A transformation extending the result 
in symmetric direct graphs to undirected graphs is established in section V. In this way, upper bounds of the number 
of Hamiltonian cycles in undirected graphs is also obtained. Some concluding remarks are proposed in section VI. 



II. NATURE BOUNDS VIA MATRIX PERMANENT 



To establish the results on the number of Hamiltonian cycles, some related concepts and results are introduced. 
Consider G — (V,E) be a directed graph with vertex set V = {vi,V2, ■ ■ ■ ,v n } and edge set E. In the following 
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section, our notations are only related to the directed graph except in section V. 

Definition 2.1 Directed graph G is called a symmetric directed graph iff edge (Vi,Vj) G E =>• (vj,Vi) G E, 
for any i =/= j. 

Definition 2.2 An 1-factor of a directed graph G is a spanning subgraph of G in which all in-degrees and 
out-degrees are 1. 

An example of the 1-factor is a spanning union of vertex disjoint directed cycles. Let NH(G) and F(G) denote the 
number of Hamiltonian cycles and the number of 1-factors of a graph G respectively. Since every Hamiltonian cycle 
is also an 1-factor, therefore 

NH(G) < F(G) (1) 
The permanent of a matrix A — (aij) nxn is defined as 

n 

Per(A) =J2H a M>) ( 2 ) 

<T 1=1 

where a goes over all the permutations {1, 2, • • ■ , n}. The adjacency matrix A = Aq of a graph G is an n by n 0-1 
matrix. The matrix A = (a^) is defined as a^j = 1, if (vi,Vj) G E; = 0, otherwise. Note the diagonal entries of A 
are all zero, and for any permutation a, a ia ^ = 1, i = 1, • • • , n iff their corresponding edges in G form an 1-factor of 
G. Hence, 

F(G) = Per(A). (3) 

Hence any upper bounds on matrix permanent would provide upper bounds for the number of 1-factors and therefore 
the number of Hamiltonian cycles. For the permanent of a matrix, the following are two famous upper bounds. 



We present a similar definition as permanent of a matrix, and we called it Hamilton of a matrix, which is 
defined as ham(A) = an, when n = 1, and when n > 2, 

ham(A) = ak 1 k 2 a k 2 k 3 ■ ■ • afc„_ifc„afe„fei , 

{kz,k3,- - ,kn} 

where {k2,ks,- • • ,/c n } is over all the permutations of {1,2,--- ,n}/{fci} and k\ is any number from the set 
{l,2,---,n}. 

Considering the relation between A and graph G, the elements in the set {ak 1 k 2 7 a k 2 k 3 , ■ • •afc 7l _ 1 fc n ,a / t 7i fc 1 } are all 
positive iff their corresponding edges in G form a Hamiltonian cycle. Hence, ham(A) = NH(G). 

Theorem 2.1 (Mine's Bound) Let A — (a,i.j) be an n x n 0-1 matrix with the row sum i — 1,2, ••• ,n. 
Then 

n , 

Per(A)<l[ r -^. (4) 

i=l 

Theorem 2.2 (Bregman's Bound) Let A = (aij) nxn be an n x n 0-1 matrix, and 7"j denote the number of ones in 
the row i, 1 < z < n. Then 

n 

Per(A) < (5) 
»=i 

Due to Stirling formula n\ < V 27rnn"/e™~ 1 / 12 , the bound in Theorem 2.2 is tighter than that in Theorem 2.1. The 
bound in Theorem 2.2 is conjectured by Mine in 1963 [14| and later proved by Bregman[5]. It plays an essential role 
in the proof of the conjecture of Szele by N. AlonQ. From the formula (1)(3)(4)(5), we can naturally obtain an upper 
bound of the number of Hamiltonian cycles of a directed graph. Particularly it is an upper bound of the number of 
Hamiltonian cycles of a symmetric directed graph. 
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III. MODIFIED RASSMUSSEN'S ALGORITHM 

In this section, We suppose G = (V, E) be a directed graph with n vertices {1, 2, • • • , n} and A = Aq = (a,ij) n xn 
be the adjacent matrix of the graph G. Let A(ij) be the (n — 1) x (n — 1) matrix obtained by removing row i and 
column j from the matrix A; A(i, :) is row i of the matrix A. For any set S, let |5| be the number of its elements. 
We now present the algorithm given by Rassmussen in [l5j . 

Algorithm 3.1 

inputs: A, an n x n 0-1 matrix; 

outputs: Xa, the estimator of the number of Hamiltonian cycles in G; 
stepO: Let pi = 0, i = 1, • • • , n; 
stepl: For i = 1 to n 

If \A(1, :)| = 1; Set p n = an, goto step2; 
Else W = {j > 1 : ay = 1}; 
If W = 0; Set A A = 0; 

Stop; 
Else 

Choose J from W uniformly at random; 
Let Pi = \W\; 

Permutate the column 1 and J; 
Let A = A(ll); 
step2: Xa = p% k ■ ■ ■ x p n . 

Algorithm 3.1 presents an unbiased estimator of the number of Hamiltonian cycles G, which means the expectation 
of the output X a is the number of Hamiltonian cycles in G. We present another point of view of Algorithm 3.1. 
Through one random experiment of Algorithm 3.1, if the output Ya is not zero, one obtains a Hamiltonian cycle of 
G by putting together all the edges corresponding to the selected elements in A. Hence, each Hamiltonian cycle of G 
can be selected with certain probability. Suppose the Hamiltonian cycle (1, fci, • • • , fc n _i, 1) has been chosen, where 
{fcj, i = 1, 2, • • • , n — 1} is a permutation of {2, 3, • • • , n}. In ith iteration of stepl in Algorithm 3.1, the probability 
of some element a%j is selected with probability 

— , i = l,--- ,n. 

Vi 

Hence this corresponding Hamiltonian cycle formed by the edges corresponding to the chosen elements is selected 
with the probability — j- x ■ • • x — and the output is p\ x ■ • • x p n . From this viewpoint, the output is an unbiased 
estimator of the number of Hamiltonian cycles of G. We state it below. 

Theorem 3.1 Let X A the output of Algorithm 3.1 . Then E{X A ) = NH(G). 

Proof: Let 7i(i) be one selected Hamiltonian cycle and denote the output when H(i) is selected. From the 

above analysis, we see Tt(i) can be selected with probability , Hence, 

NH(G) 

E(X A )= J2 Y H{l) — = NH{G) 

gives the result. □ 

Note that in one random experiment in Algorithm 3.1 in stepl, the element is selected by the ascending order of 
the row, or equivalently say, we select the first element from the first row of A, then select the second element from 
the second row of A and so on. If we select the element by another fixed order of the row, we can obtain Algorithm 
3.2, which performs an essential part in getting new upper bound in a symmetric directed graph. In the following 
Algorithm 3.2, the matrix B is used to determine which row are selected at each iteration stepl of Algorithm 3.2 in a 
random experiment to construct a Hamiltonian cycle. In each independently random experiment running Algorithm 
3.2, if B is chosen, it remains unchanged, which promises the results from the independent experiment to be identical 
random variables. 

Algorithm 3.2 
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inputs: A, an adjacent matrix of graph G; 

B = (bij) nxn a matrix, where 6^ is chosen from {1, 2, • • • , n — i + 1}, for any i,j € n. 
outputs: Xa, the estimator of the number of Hamiltonian cycles in G; 
stepO: Let pt = 0, i = 1, • • • , n and fc = 1. 
stepl: For i = 1 to n 
Set gi = bik ; 
If 1^,01 = 1; 

Set p n = a 3i i, goto step2; 
Else W = {j + 9l : a gd = 1} ' 
If W = 0; Set X A = 0; Stop; 
Else 

Choose J from W uniformly at random; 
Let pi = \W\ and fc — J; 
Permutate the column gi and J; 
Let A = A(gigi); 
step2: Xa = Pi x • • • x p„. 

We now show that Algorithm 3.2 presents an unbiased estimator of the number of Hamiltonian cycles in G. 
Before present the proof we need a technique lemma. 

Lemma 3.1 Let A'(ij) denote the matrix obtained from A by removing row i and column i after pcrmutat- 
ing column i and column j. Then, if n > 2 

ham(A) = ^ a^j h&m(A'(kij)). 



Proof: Without loss of generality, we suppose fci = 1. 
The base case n = 2 is trivial. 
Suppose the case n — 1 holds. 
Now we see the case n. Since 



By induction, we know 



ham(A) 



^ aifc 2 a fc 2 fc 3 • • • afe rl _ 1 fe li afc re i 

n 

aij<ijk 3 ■ ■ ■ ak n - l k n a-k n \ 

j=2 {k 3 ,- ,k n } 
n 

X/ ai J a J*=3 ' ' ' a fe„-l*:„ a fe„l 

j=2 {fc 3 ,-,M 



where {fc 3 , • • • , fc„} is over {2, • • • , 
Hence, we need only to show 



a J fe 3 ' ' ' afc„_ lfc „a fc „i = ham(A'(lj)) 

{k 3 ,-,k n } 



ham(A'(lj)) = £ 



fc' fc' fc' fc' 



■ ' a k' n _ 1 k , n a k' n i 



where {k' 3 , • • • , k' n } is over all the permutations of {2, • • • ,n— 1}. Recall the definition of A'(lj), which is obtained 
by removing row 1 and column 1 after permutating column 1 and column j, then we know a' lk , = aj t k' + i and 

a' k , +1 k , +1 , a' k , 1 = a' k , which completes the proof. □ 



A. o fc ^ fc q "I - 1 t fc i "I - 1' ' fcrj 1 ^ ■ 
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By lemma 3.1, it's sufficient to show Algorithm 3.2 presents an unbiased estimator of the number of Hamiltonian 
cycles. 

Theorem 3.2 Let X A the output of Algorithm 3.1. Then E(X A ) = NH(G). 

Proof: We go on to show E{Xa) = NH(G) by induction on n. 

The base case n = 1 is trivial. 

Suppose the case n — 1 holds, for the case n, 

E(X A ) = E(X A \J = j)P(J = j) 
= E(piX A , {klj) )/ Pl 

= E( y X A'(k 1] ))- 

j¥*;i 

By induction, E(X A ^ kl j- ) ) = ham(A'(fci j)) and Lemma 3.1, then 

E(X A ) = ham(A , (fc 1 j)) = ham(A) = NH(G). 

j¥fei 

This completes the proof. □ 

IV. AN UPPER BOUND OF SYMMETRIC DIRECTED GRAPHS 

We now present the upper bounds for the symmetric directed graph. 

Theorem 4.1 Let G be a symmetric directed graph, A = Aq = (aij) n xn be the adjacent matrix of G and 
Ti denote its sum of row i, i = 1, 2, • • • , n, n > 3. N denotes the number of Hamiltonian cycles of G. Then 

n 

N < 7: — rY\n 

Proof: Let Tt(j) = (1711,1712, ■ ■ ■ , m n , mi) be one of the Hamiltonian cycles of G, where j = 1, 2, • • • , N. mi (i = 
2, 3, • • • , n) is a permutation of {1,2,- • • ,n}/{mi}. In Algorithm 3.2, choose b\\ — m\. In step 1, choose elements 
by the row order of m 2 , - • • , m n such that the edges corresponding to the chosen elements constituting H(j). Let 
Si = {j : a mi j — 1} and "pf denote the value of ith iteration of stepl in Algorithm 3.2, where i = 1,- • • ,n and 
m n+ i = mi. Then 

Pi = r mi j 

vf = \Si/{m 1 ,m 2 , ■ ■ ■ i = 2, • • • ,n- 1, 

^ = 1- 

n 

Let Xft(j) be the output when the edges corresponding to the chosen elements form H(j) , then X n ^ = \\ pf . 

i=l 

Since this is a symmetric directed graph and n > 3, there exits a Hamiltonian cycle 7i (j) = (mi, m„, m n _i, • • • , mi) 
different from H(j). Let pf be the value of ith iteration of stepl in Algorithm 3.2 and X n ,^ the output when the 
edges corresponding to the chosen elements form H (i), where i = 1, • • • ,n and m = m n . Then by Algorithm 3.2, 

Pi = f mi , p§ = 1 , 

pf = \Si/{mi,mi,m i+1 , ■ ■ ■ ,m n }\, i = 3,--- ,n, 

n 

*wu) = n pf ■ 

Therefore 

n 

X H{j)X H '{j) = \\pfpf ■ 
i=l 



6 



Considering the symmetry of A and a mimi = 0, i = 1, 2, • • • , n, then pf+pf < r mi , i = 2, • • • , n. Hence 

n 

x nU) x H , U) = Y[ P ?pf 

i=l 

n 

<r 2 mi Hpf(r mi -vf) 

i=2 
1 ™ 

— TTr 2 

n-1 11 rm 



< 



4_n- 



i=l 



n 

-—Y[r 2 

4n-l 11 ' 
i=l 



Let X A be the output of Algorithm 3.2. Then by Theorem 3.1 and Theorem 3.2 P(X A = X n t^) 



From 



N 



Algorithm 3.2, we know the output may be zero with certain probability, hence ^ x 1 < 1 . Set N — NH(G), we 



have 



i=i 



HO) 



N < 



< 



N 



N 



N 



n x w) 



N 



\ H x n(j)X W (j) 



< 



IV IL IL 

^ n 4n^i n r ? = n r *- 



Thus the result follows. □ 



Theorem 4.2 Let A = (a,ij) nxn be an adjacent matrix of a symmetric directed graph and be the sum of 
row z of A. Then 



_^ n n 

— n^n 



r, + 1 



i=i i=i 



Proof: Since 



2=1 



>I(1 + I r >l. 



Thus the result follows. □ 



Theorem 4.2 shows our upper bound is tighter than Mine's bound (4). In many cases, the new bound is better 
than Bregman's bound (5). For example, 



n n 



i=i i=i 



when n > 100, r» < 5, i = 1, 2, • • • , n. 
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V. BOUNDS OF UNDIRECTED GRAPHS 

The notations or definitions related to the undirected graphs arc only stated in this section. The problem of 
counting the number of Hamiltonian cycles in an undirected graph is transformed to that of counting the number of 
Hamiltonian cycles in a symmetric directed graph. This transformation is very simple but powerful. Let G be an 
undirected graph with vertices {1, 2, • • • , n}, where n > 3. G is a simple graph. Define a symmetric directed graph 
G' corresponding to G by replacing each edge of G with two directed edges and Let Hq and Hqi 

denote the set of the Hamiltonian cycles in G and G' respectively. V{Hqi) denotes the power set of He- Recall we 
use the notation v\v 2 ■ ■ ■ v n v\ and (vi,V2, ••• ,v n , v\) to denote Hamiltonian cycles in undirected and directed graphs 
respectively. 

Theorem 5.1 Let H = m\m<2, ■ ■ ■ m n m\ be a Hamiltonian cycle in Hg- Then there are at least two Hamil- 
tonian cycles (mi,TO2,-" - , m„,mi) and (mi, m„, m n _i, • • • ,mi) in He- Define a map ip from Hq to V(Hc>) as 
follows: 

<p(H) = {(mi,m 2 , • • • ,m n ,mi), (mi,m„,m n _i, • • • ,mi)}. 

Let Imtp denote the image of the map if and H' = m' x m'^ ■ ■ ■ m' n m\ be a different Hamiltonian cycles from H. Then 

ip(H) ("1 ip(H') = and {Jlrmp = H G >. 

Proof: Due to the symmetry of the graph and n > 3, if there is a Hamiltonian cycle (mi,m2, • • • , m„,mi) in i?G') 
there must be a different Hamiltonian cycle (mi,m n , m„_i, • • • , mi) in These two Hamiltonian cycles obviously 
has a pre-imagine corresponding to the Hamiltonian cycle m,\m,2 • • • m n m\ in Hq- Note (mi, m,2, ••• , m„, mi) is in 
ip(m\m2 ■ ■ ■m n mi). Hence, Ulnvp 2 Obviously, Ulrmp C iJc- Therefore 

Lilmtp = He ■ 

Suppose there are two different Hamiltonian cycles W = m\m2 ■ ■ ■ m n m\ and W = mjm^ ■ • • m! n m! 1 in We 
know that they are different iff there exits a vertex {m,}={m'} such that at least one of the two neighbor 
vertices of {m^} is not in the set of two neighbor vertices of {m 1 -}. Hence (mi,m2, • • • ,m n ,m{) is different from 
(m'i,m 2 ,--- ,m^,m'i) and (m' 1 ,m' n ,m' n _ 1 ■ ■ ■ , m^mi), we know (mi,m 2 , ■ ■ ■ , m„,mi) is not in the set <p(H r ). 
Similarly, (mi, m„, m„_i, • • • ,mi) is not in f(H'). Then <^(H) n ¥?(W) =0. □ 

Corollary 5.1 Let NH(G) and NH(G') denote the the number of Hamiltonian cycles in undirected graph 
G and its corresponding symmetric directed graph G' respectively, then 

NH{G) = ^NH(G') 
Proof: This result is a straightforward deduction of the Theorem 5.1. □ 

Corollary 5.2 Let G be an undirected graph with vertices {1,2,- •• , n}, n > 3, NH(G) be the number of 
Hamiltonian cycles in the graph G. di denotes the degree of the vertex {«}. Then 

1 ™ 

z i=i 

n 

A^(G)<-n^ 

^ i= i 

and 

1 ™ 

^(G)<-n(^) iM - 

i=l 

Proof: By the results of Theorem 2.1, 2.2, 4.1 and Corollary 5.1, this corollary follows. □ 
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VI. CONCLUDING REMARKS 



A novel upper bound of the number of Hamiltonian cycles on a symmetric directed graph is presented first in this 
paper, which is tighter than the famous Mine's bound and better than the bound by Bregmman in many cases. 

A transformation from the problem of Hamiltonian cycles of an undirected graph to that of the symmetric directed 
graph is constructed. Using this transformation and the bounds for directed graphs, upper bounds for the number of 
Hamiltonian cycles in undirected graph are obtained. The significance of this transformation also lies in the fact that 
the algorithms for counting the number of Hamiltonian cycles in a directed graph can be directly applied to count 
the number of Hamiltonian cycles in an undirected graph. 
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